Applications of four-body exponentially correlated functions 
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We demonstrate the applicability of four-body exponentially correlated functions for the accurate calculations 
of relativistic effects in lithium-like atoms and present results for matrix elements of various operators which 
involve negative powers of interparticle distances. 
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I. INTRODUCTION 

In the accurate evaluation of atomic energy levels not only 
the nonrelativistic energy, but also relativistic and QED ef- 
fects have to be calculated with the high precision. The only 
approach which consistently accounts for all corrections in 
small atomic systems is the one based on the expansion of 
the energy in the fine structure constant a 
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m a 



(4) 
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Each term in this expansion can be expressed as the expec- 
tation value of some effective Hamiltonian with the nonrela- 
tivistic wave function. Namely f^'^' is the nonrelativistic en- 
ergy (in atomic units), is the relativistic correction, which 
for states with the vanishing angular momentum is the expec- 
tation value of i/*^'*' in Eq. f'^^ and higher order cor- 
rections are expressed in terms of matrix elements of some 
more complicated operators. They have been calculated for 
low lying states of the helium atom and helium-like ions up to 
the order ma^ [1], and for the particularly important case of 
2^Pj splitting up to the order m f^]. One of the sources of 
this achievement was the flexibility of the explicitly correlated 
exponential basis set, which due to its correct analytic prop- 
erties, makes possible accurate evaluation of matrix elements 
with complicated and singular operators. 

In the case of the lithium atom and light lithium-like ions 
all corrections up to f have been accurately calculated fsL 
m, but not that of higher orders. The principal reason for the 
much slower progress for three-electron atoms is the difficulty 
in handling integrals with explicitly coiTelated functions. The 
commonly used explicitly correlated Gaussian functions do 
not have right analytic properties, for example they do not 
satisfy the cusp condition, and therefore cannot be used for the 
calculation of higher order relativistic corrections, like 
Hylleraas basis functions have the right analytic behavior: the 
accuracy in solving the Schrodinger equation is the highest 
among all other basis functions, but it is difficult to handle 
Hylleraas integrals involving quadratic negative powers of two 
different interparticle distances. Such integrals appear in the 
evaluation of f and for this reason other basis functions 
have been investigated in the literature. 



■Impuchals @ fuw.edu.pl| 
fkrp@fuw.edu.pl| 



Zotev and Rebane in were the first to apply exponen- 
tially correlated functions 
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in variational calculations for Ps2 and the other exotic 
molecules. They have found a simplified formula for ma- 
trix elements of the nonrelativistic Hamiltonian and presented 
numerical results of variational calculations with a few basis 
functions. In our recent paper |7] we presented an efficient 
algorithm for the evaluation of integrals involving powers of 
ri and r;,- 
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with Ua being nonnegative integers. It is based on recur- 
sion relations which start from the master Fromm-Hill integral 
lUllst], where all Ua = 0. We have applied this algorithm to 
the variational calculations of the ground state of Li and Be+ 
with up to 128 functions. The comparison of nonrelativistic 
energies with the ones obtained with much larger number of 
Hylleraas functions indicates that the exponential representa- 
tion of the three-electron wave function is very efficient. 

The class of integrals in Eq. ^ with nonnegative Ua is 
sufficient for nonrelativistic energies (0, ^ . However, it is not 
sufficient to calculate the leading relativistic effects described 
by Breit-Pauli Hamiltonian, which for S-states takes the form 
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Its matrix elements involve an extended class of integrals with 
exactly one of ria equal to —1, and all others are nonnegative, 
while that for leading QED effects involve integrals with Ua — 
—2. This is well known from calculations with Hylleraas basis 
functions, where all Ua in Eq. Q are equal to zero. Hylleraas 
extended integrals of that kind have been extensively studied 
in iflol - ilTll using multipole-type of expansions and rec ent l y b y 
present authors using analytical recurrence approach 
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Both methods had been successfully applied in high-precision 
calculations of leading relativistic and QED corrections to the 
energy of lithiumlike systems ll3[ [l9ll20t] . There are no similar 
studies for exponentially correlated integrals to the best our 
knowledge, and for the first time we present them in this work. 

In the calculation of relativistic and QED effects beyond 
leading order, f^^^ for example, another class of integrals ap- 
pears with two quadratic inverse powers of interparticle dis- 
tances. There are only few studies in the hterature for three- 
electron Hylleraas integrals lfl2l [Tsi [TtIi . The algorithm by 
King [ 17] seems to be too slow for a large scale computation, 
where integrals with il = J^a "a °f order 30 have to be per- 
formed. The evaluation of these integrals is quite difficult with 
the recursion method and have not yet been worked out so far. 

In the case of exponentially correlated integrals the prob- 
lem seems to be even more severe, since the master integral 
with Ua ^ is much more complicated. However, being 
able to optimize nonlinear parameters of each function inde- 
pendly, one does not need to use large powers of interparti- 
cle distances in the basis set. For S-states it is sufficient to 
use functions of the form (|2]i. In such a case, having an ana- 
lytical and thus accurate method for g{ni, n2, ns, n4, ns, uq) 
with nonnegative ria, inverse negative powers of the interparti- 
cle distance can be obtained by the numerical integration with 
respect to the corresponding parameter Wa or Ua- It however 
requires a good control of numerical accuracy of the master 
integral and of recursion relations in Eq. (fTOl i. The usage the 
higher precision arithmetic is essential in some critical areas 
of the integration. A key feature of our numerical integration 
strategy is the adapted quadrature, which allows one to get the 
high accuracy with a very small number of points. 

We demonstrate our method on examples with expectation 
value of various operators on lithium ground state. Results 
obtained for matrix elements involving single Ua = —1 are 
compared to the most accurate ones obtained with the Hyller- 
aas basis set. Good numerical convergence of results for 
matrix elements involving two negative powers, for example 
^/(rlrl), ^/{rlrl^), l/irl^rlJ. indicate that this integration 
approach can be used for the calculation of higher order rela- 
tivistic corrections, for example m and m effects in the 
hyperfine and fine structure of lithium-like systems. 



group 5*4. We can assign vertices 1,2,3 to the electrons and 
to the nucleus as shown in Fig. ([T), and to edges we assign 
Ua,Wa and Ua parameters of a given integral. The symme- 
try group 54 corresponds to 24 renumbering of vertices 0,1,2, 
and 3, which means also a relevant change of parameters on 
the edges. The generated symmetry relations allow us to re- 
duce the number of recurrence formulas for Slater integrals. 
It is necessary to derive only one recurrence scheme, and the 
formulas for the advancement in the other indices can then be 
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FIG. 1 : Tetrahedron representing the integral from Eq. (|2} 

obtained by application of the ^4 symmetry. 



B. Basic integrals 

The evaluation method for g{ni, 712, n^, 714, n^, ng) for the 
basis class with all nonnegative was first presented by Har- 
ris in [21] and later by us in Ref. [7]. Here we present only a 
short summary, which is needed for the evaluation of the ex- 
tended integrals. The master integral go = g{0, 0, 0, 0, 0, 0) 
satisfies the following differential equation 
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where the 5*4 symmetric polynomial a is of the form 
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II. CALCULATION OF INTEGRALS 

A. Tetraliedral symmetry 

An important property of integrals defined in Eq. (|3]l is the 
tetrahedral symmetry which is equivalent to the permutation 
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P = -UiWi[{ui + W2)'^ - ul]r{0,0,-l;Ui+ W2,U3,U2 + Wi) 

-ui wi [{ui + U'if - wl] r(0, 0, —1; Ml + -U3, 102, wi + 103) 

+ [ul wl -\-u\w2 ~ u\ w\ + Wl W2 {u\ +u\— w1)\ r(0, 0, -1; Wi + W2, W3, Ul + U2) 
-\-[u\ w\ - u\w2+ u\ wl + Wl W3 {uf + U3 — wl)] r(0, 0, -1; Wi + W3, W2,Ul + U3) 

-[u2 {u2 + Wl) {u\ + U3 — w\) - u\ {ul + u\~ wl)] r(0, 0, -1; U2 + "011,^3, ui + ^2) 
-["3 ("3 + Wl) {u\ +ul- wl) - ul {u\ + U3 - ^2)] r(0,0, -l;u3 + wi,U2,ui + W3) 
+wi [w2 {u\ -ul+ wl) + W3 {u\ +wl- ul)] r(0, 0, -1; W2 + ^3, Wi,M2 + -"3) 

+Wi [U2 [ul -~wl+ ul) + U3 [ul +ul~ wl)] r(0, 0, -1; U2 + U3, Wi,W2+ W3) , (7) 



and where the two-electron integral F is 



r(ni,n2,n3,a,^,7) = 
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The recurrence relations for the integrals g result from the dif- 
ferential equation (|5]l and can be written in a compact form 
as 
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The relation ( fTOl i allows one to express the integral 
g(ni, ..jrig) with positive index ni through g-integrals with 
smaller nonnegative indices. Derivatives of polynomials 
in P(ni,n2, "-3, "-4,^5,716) and cr(ni, 712, 713, 714, ns, ng) are 
calculated explicitly. The master integral g(0, 0, 0, 0, 0, 0) and 
the two-electron integrals F are needed to start the evalua- 
tion of g(77i, 7i2, 713, 714, 715, 7i6) from recursion relations. For 
0(0, 0, 0, 0, 0, 0) we implemented the formula of Fromm-Hill 
dSD in the version improved by Harris [9]. The calculation of 
the two-electron F functions has been described in detail in 
Refs. I,22r,24 1 . In Eq. (fTOl i the parameter wi is distinguished 
on the right-hand side. We can define the same recurrence re- 
lations with other variables Wa , Ua from expressions (Q, ( fTOl l. 
and (fTsT i by applying the tetrahedral symmetry. 



where we use a Newton-like symbol notation 
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P(ni, 712, 713, 714, rig, Tig) = 
(T(7ll,7i2,"-3,"-4,"-5,"-6) 
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The proposed recurrence scheme allows us to calculate in- 
tegrals from higher shells il very efficiently. In Table |T] we 
present values in two reference points introduced by Fromm 
and Hill in Ref. fl. These are the standard reference point 
(SRP) Wa — Ua = 1 where <t = —2 and the auxiliary refer- 
ence point (ARP) Wa — l,Ua ~ with <t — 1. Values for the 
last one can be compared to the known Hylleraas results. All 
presented digits are significant, which confirms the very good 
stability of the recursions at least at these reference points. 



The integrals g{ni,n2, 713, 714, 715, Tig) with cr is close to zero 
are difficult to evaluate because the recurrence relations are 
numerically unstable. Following Harris' studies on the mas- 



ter integral g{0, 0, 0, 0, 0, 0) in Ref. we considered Wa — 
l,Ua — a close to a = 3^^/^ where a is exactly equal 
to zero. As an example we present in Table |II| values for 
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TABLE I: Exponentially correlated integrals in SRP and ARP reference points 



ARP : tia = 0, Mia = 1 



SRP : M„ = to„ = 1 



g(0, 0,0, 0,0,0) 
0,0, 0,0,0) 
1,0,0,0,0) 
9(2,0,0,0,0,0) 
g(l, 1, 1,0,0,0) 
9(3,0,0,0,0,0) 
9(1,1,1,1,0,0) 
9(4,0,0,0,0,0) 
9(1,1,1,1,1,0) 
9(5,0,0,0,0,0) 
9(1, 1, 1, 1, 1, 1) 
9(6,0,0,0,0,0) 
9(2, 2, 2, 2, 2, 2) 
9(12,0,0,0,0,0) 
9(3,3,3,3,3,3) 
9(18,0,0,0,0,0) 



2.208 310 154 388 618 874 536 424 144 10"^ 
2.208 310 154 388 618 874 536 424 144 10"^ 
2.587 328 950 655 650 145 088 210 000 9 10"^ 
3.658 582 716 243 175 207 969 276 574 1 10"^ 
3.447454 259 102 525 282 100 115 713 9 10"^ 
8.803 723 087 040 150 596 505 448 579 7 10"^ 
8.518 518 518 518 518 518 518 518 518 5 lO" 
2.849 464173 126 685 211 199 798 289 3 10° 
2.500 000 000 000 000 000 000 000 000 10° 
1.176 795 411 671 425 935 279 581 659 10^ 
8.000 000 000 000 000 000 000 000 000 10° 
5.962 899 567 501 152 778 486 008 087 6 10^ 
4.761 568 658 224 658 953 775 935 886 3 10* 
1.998 667 989 835 473 493 169 581 298 3 lO'' 
4.595 097 600 000 000 000 000 000 000 10° 
1.779 533 879 212 729 823 179 024 201 1 10^* 



2.173 757 633 275 034 284 081 325 213 9 10" 
1.086 878 816 637 517 142 040 662 607 10" 
7.126 296 734 072 723 436 203 022 296 410" 
8.571 078 153 650 184 467 645 704 204 7 10" 
6.027 278 683 541 261 617 671 441 399 9 10" 
9.316 428 522 373 859 205 286 070 304 7 10" 
4.807 529 762 877 014 687 398 150 589 3 10" 
1.296 572 573 872 689 431 790 047 714 4 10" 
4.679 985 481 131 151 729 598 237 602 7 10" 
2.205 496 365 095 795 104 480 271 164 6 10" 
5.033 124 034 041 977 640 019 585 799 2 10" 
4.441 445 631 268 318 817 661 363 268 6 10" 
1.467 676 925 382 090 748 693 886 720 7 10" 
3.793 674 032 518 970 385 528 210 586 4 10^ 
5.979 770 415 024 714 287 528 834 619 9 10^ 
6.240 107 562 655 131 418 896 693 018 6 10^ 



TABLE 11: Numerical values of g{l, 1, 1, 1, 1, 1) around the singular point a = 0: 101,2.3 = 1, wi 



2,3 ~ OL- 



a 



1.741 3 10"* 
2.396 798 8 10"^ 
6.568 82917 10"^° 
8.925 087 775 552 10" 
2.648 337 377 078 125 10" 
3.169 230 531337161 925 10" 
2.788 669 608 438 596 620 038 649 3 10" 


-6.754 320 066 991579 670 162 432 10" 
-3.147178 562 004 038 394 8 10"^° 
-8.157642 380 596 28 10" 
-2.571592 837 582 7 10" 
-2.807218 7 10~° 
-1.067302 7 10"^ 
-1.722 8 10"'' 



9(1,1,1,1,1,1) 



4.715 168 989 550 199 879 343 821 168 10"^ 
4.713 684 016 763 871 477 743 039 334 10"^ 
4.713 681 976 034 959 469 113 297 912 10""^ 
4.713 681 970 427 392 724 217 475 890 10"^ 
4.713 681 970 426 653 329 868 706 104 10"^ 
4.713 681 970 426 630 719 189 520 726 10"^ 
4.713 681 970 426 630 718 918 940 842 10"^ 
4.713 681 970 426 630 718 918 938 462 10"^ 
4.713 681 970 426 630 718 918 937 885 10"^ 
4.713 681 970 426 630 716 231 943 331 10"^ 
4.713 681 970 426 623 754 094 755 315 10"^ 
4.713 681 970 424 435 146 822 397 542 10 
4.713 681 946 459 185 584 378 695 888; 10 
4.713 681059 186 068 003 930 086 238 10"^ 
4.712211406 354 297107 821 210234 10"^ 



0.5773 

0.577 350 2 

0.577 350 269 

0.577 350 269 189 6 

0.577 350 269189 625 

0.577 350 269 189 625 764 5 

0.577 350 269 189 625 764 509 148 7 

3.-1/' 

0.577 350 269 189 625 764 509 148 8 

0.577 350 269 189 625 764 6 

0.577 350 269189 626 

0.577 350 269189 7 

0.577 350 270 

0.577 350 3 

0.5774 



-,-15 



-,-20 



,-16 
-13 



-2 
-2 



(7(1,1,1,1,1,1) which needs six evaluations of recursions 
( [Tol l. Close to the critical point a = 0, we used Bailey's mul- 
tiprecision library |25]. We control the number of significant 
digits by dynamical estimating and adjusting the appropriate 
precision of the arithmetic. As presented in Table [III we can 
approach the critical point cr = as close as we need for prac- 
tical purposes. This strategy of course slows down the algo- 
rithm significantly, but in practical applications the parameters 
close to the critical point = are very rare. This strategy of 
controlling precision in the region of instabilities allows one 
to cross a — points in the minimization of the nonrelativis- 
tic energy. In Table |III] we present results for nonrelativistic 
energies for the ground state of Li obtained with global min- 
imization of all nonlinear parameters in the basis length of 
iV = 128, 256, 512 respectively. The achieved precision is 
much higher than that from similar number of Hylleraas func- 
tions. 



TABLE 111: Nonrelativistic energy of the Li ground state using expo- 
nentially correlated basis set with the number of functions A*', oo is 
the extrapolated value from Ref. |5] using Hylleraas basis functions. 



N 




32 


-7.478 059 40 


64 


-7.478 060 050 


128 


-7.478 060 304 6 


256 


-7.478 060 32131 


512 


-7.478 060 323 427 


oo 


-7.478 060 323 9101(3) 



opportunity to obtain extended integrals by using either a sin- 
gle integration over Wa or Ua i.e. 



C. Extended integrals 

In this section we present an algorithm for calculations of 
extended integrals with 1/r^ or l/r^^ factors in Eq ([3]l. This 
means that some of indices in g(ni, n2, na, n4, ns, ng) are 
equal to ~1. Fully correlated exponent in Eq. ^ gives the 



g{ni,n2,n3, -1,715, Uq) 

/ dwiy(ni,n2,n3,0,n5,n6) , 

(14) 
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TABLE IV: Examples of extended integrals with 1 / 7-23 in SRP and ARP reference points calculated using numerical integration with A'^ = 30. 







A D . > . n .... 1 
Aixr . Ua — U, Wa. — 1 


Oixr^ : Ua — Wa — J- 




g{0, 0, 0, - 


■1,0,0) 


3.852 610 933 969 379 240 110 048 369 9 10^^ 


9.496 501 144 947 432 180 784 237 909 6 10~ 


'2 


g(l, 0, 0, - 


■1, 0, 0) 


3.Uz7 44y lOb o7d 050 3d7 30o YOz 234 4 10 


o Tor\riTo i£?r\Tcr^ /ioTr\co i>r\o or^o r\oi T ir\ — 

3.729 972 160 750 497 052 308 263 021 7 10 


2 


s(l, 1,0, - 


■1, 0, 0) 


2.360 700 062 552 231 395 696 277 796 4 10"^ 


1.895 953 170 412 559 215 217 374 496 4 10" 


2 


9(2, 0, 0, - 


■1, 0, 0) 


4.360 947 194 620 688 310 533 551 110 5 10~^ 


2.519 224 983 963 737 154 014 628 472 2 10" 


2 


s(i, 1, 1, - 


■1, 0, 0) 


2.503 315 976 630 860 793 316 473 181 3 10~^ 


1.453 670 463 091 034 929 350 269 890 4 10" 


2 


S(3,0,0, - 


■1,0,0) 


9.678 534 014 755 149 691 745 602 065 7 10"^ 


2.452 489 069 429 084 805 380 219 451 1 10" 


■2 


S(l, 1, 1, - 


■1, 1,0) 


4.965 970 761 362 395 801 363 888 667 7 10"^ 


8.929 779 730 282 261 688 941 298 690 9 10" 


3 


S(4,0,0,- 


■1,0,0) 


2.994 734 857 888 708 511 724 418 787 7 10° 


3.140 438 234 282 031 347 870 179 278 8 10" 


2 


S(l, 1, 1, - 


■1,1,1) 


1.333 333 333 333 333 333 333 333 333 3 10° 


7.833 365 875 543 837 545 344 094 254 1 10" 


3 


g(5,0,0, - 


■1,0,0) 


1.207 413 986 734 661 793 622 502 357 9 10^ 


5.003 447 967 255 362 321 317 326 500 6 10" 


2 


S(2, 2, 2, - 


■1, 2, 2) 


5.436 536 048 634 697 021 325 813 246 7 10^ 


6.776 488 986 538 662 624 450 203 053 8 10" 


2 


g(10,0,0, 


-1,0,0) 


1.820 037 338 296 110 774 925 172 748 9 lO'^ 


4.998 826 524 680 257 565 021 055 520 2 10° 




3(3,3,3,- 


■1, 3, 3) 


2.636 015 376 623 376 623 376 623 376 6 lO'' 


6.615 220 449 456 842 761 016 682 847 5 10° 




g(15,0,0, 


-1,0,0) 


4.364 311 343 749 328 364 240 301 511 7 10^° 


6.183 347 201 330 971 024 067 326 468 8 10^ 





or a double integration i.e. 

g(ni, -1,713, -1,^5,716) = 

/ dui dw2 g{ni,0,n3, 0,715, ne). 

(15) 

The adaptive increase of the arithmetic precision close to crit- 
ical points (7 = is necessary here for the precise evalu- 
ation of 712, n3, 71,4, ns, ng)- Moreover, one is able to 
perform accurately this integration by using A^-point general- 
ized Gaussian quadrature with logarithmic end-point singular- 
ity [26|] 

I = f dx[Wi{x) +\n{x)W2{x)] 
Jo 

N 

« ^ [Wi{x,) + \n{xi) W2{x,)] , (16) 

i=l 



where Wi are regular functions on the interval (0, 1). This 
quadrature becomes exact for Wi being polynomials of max- 
imal degree — 1, and the example of 30 point quadrature, 
which is used through out this paper for calculation of mean 
values, is presented in Appendix A. 

If integration variable u in Eq. ( fT4] i is mapped into the in- 
terval (0, 1) by the following change of variable 



r du f{u) = f'dx^f(--l) (17) 
Jo ^ \^ J 



then the quadrature with logarithmic end-point singularity is 
sufficient for one-dimensional integrals in Eq. ( fT4l i. where 
= 30 quadrature allows one to obtain about 30 significant 
digits as shown in Table II VI 



The presented values for integrals are obtained at ARP and 
SRP points. Some of them can be found in the literature, for 
example g(2,2,2,-l,2,2)at ARP point (nlH, which cor- 
responds to a Hylleraas type of integral. Perfect agreement 
with those results, demonstrates high accuracy is achieved for 
the extended integrals. The proposed evaluation method fully 
relies on properties of the recurrence algorithm for basis inte- 
grals, which must be very stable on the integration path over 
the corresponding parameter. 

In comparison to the one-dimensional integral in Eq. ( fT4] i. 
the convergence of two-dimensional integral Eq. ( fTSl l with re- 
spect to the number of integration points is much worse. For 
this reason we use a slightly different mapping into (0, 1) in- 
tervals, which is 

f-OO f-OO 

/ du \ dwf{u,w) = 
Jo Jo 

[\x r dy^f(^-l,l-l). (18) 

Jo Jo x^y \^ y- J 



The numerical convergence of the integral in Eq. ( fTSl ) is the 
worst for the case = 0, where the leading asymptotics in- 
cludes a square of the logarithm. For 773 > convergence im- 
proves significantly. The use of Gaussian quadrature adapted 
to logarithmic end-point singularity with 30 points is enough 
for practical applications. In Table [Vl we presented numerical 
values for 5(771, —1, 773, —1, 715, 775) in ARP and SRP refer- 
ence points with the accuracy of 10^^^. In the case of 773 = 
they have been obtained with 60 point quadrature. It is pos- 
sible to obtain even higher accuracy for 773 > 0, but further 
improvement for 713 = requires a more sophisticated inte- 
gration strategy. There are no such problems with integration 
involving parameters which are attached to opposite edges of 
the tetrahedron, i.e. g{~l, 772, 713, —1, 715, 775), see Fig. [T] so 
this case of integral in Eq. ( fT4l i with 773 = is the one which 
limits accuracy of mean values. 
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TABLE VI: Expectation value for Breit-Pauli operators for the ground state involving integrals with one rn — —1. Implicit summation over 
a, or over pairs a > b is assumed. Last entries are extrapolated results obtained in Hylleraas basis set 



N 




n h 


"5" ('■a) 




Pa 


Pa rJiS"rl^+rl,ri^)pl 


1 


oO.VoZ 7y7 yob 7 


A T^nc /ICC crri/i i 

4.506 456 5u4 1 


13.764 569 952 4 


0.536 449 2U8 3 


625.582 840 9 


0.936 654 826 


2 


29.747 608 655 


4.576 769 565 1 


13.585 859 628 1 


0.561647 862 3 


611.553 996 7 


1.138 521450 


4 


30.130 068 846 4 


4.443 919 4210 


13.787803 1290 


0.543 639 175 1 


625.468 724 2 


0.914892 317 


8 


30.187 4817320 


4.421446 718 9 


13.815 726 640 6 


0.543 582 099 9 


627.407 849 4 


0.903179 423 


16 


30.241078 773 


4.381681063 4 


13.842 637 278 6 


0.544 318 958 8 


628.479 718 2 


0.871711023 


32 


30.240 966 286 8 


4.381283 593 8 


13.842 598 063 3 


0.544 325 804 3 


628.457 736 5 


0.871331224 


64 


30.240 892 554 9 


4.381232 0316 


13.842 567 782 2 


0.544 325 359 5 


628.451398 8 


0.871268 418 


128 


30.240 987 196 9 


4.381 186 5673 


13.842 617 080 2 


0.544 324 836 8 


628.450 904 8 


0.871208 043 


256 


30.240 973 605 8 


4.381 176 947 3 


13.842 611088 3 


0.544 324 684 9 


628.449 069 5 


0.871 196 220 


Hyll. 


30.240 972 72(3) 


4.381 176 64(4) 


13.842 610 86(3) 


0.544 324 632 5(7) 


628.448 985(12) 


0.871 195 62(14) 



TABLE V: Examples of extended integrals g{ni, — 1, na, — 1, ns, ng) in SRP and ARP reference points, numerical quadrature with N = 60 
points, all digits are significant 









ARP : Ua = 0, Wa = 1 


SRP : Ua — Wa = 1 




-1,0, 


-1,0,0) 


1.884 392 088158 216 10" 


8.424 892 130134 382 10"' 




-1,0, 


-1,0,0) 


9.485 660 506 739 961 10"^ 


1.797 709 314 546 008 10"^ 


5(1, 


-1,1, 


-1,0,0) 


4.284 596 512 743 028 10"^ 


4.614 367 395 977 269 10"^ 


5(2, 


-1,0, 


-1,0,0) 


1.165 927 539 416 184 10° 


8.712176 344 048 701 10"^ 


fl(l. 


-1,1, 


-1,1,0) 


6.829 291358121588 10"^ 


2.470 684 625 918 397 10"^ 


5(3, 


-1,0, 


-1,0,0) 


2.440 178 087422 623 10° 


6.789 663 816 285 613 10"^ 


5(1, 


-1,1, 


-1,1,1) 


1.467 401 100 272 340 10° 


1.868 910 909 646 633 10"^ 


5(4, 


-1,0, 


-1,0,0) 


7.400 703 871 489 098 10° 


7.375 782 376 091 605 10"^ 


5(2, 


-1,2, 


-1,2,2) 


1.325 932 535 285 045 10^ 


3.788 534 374 425 215 10"^ 


5(7, 


-1,0, 


-1,0,0) 


8.898 402 193 692 669 10^ 


3.628 968 910 010 360 10"^ 



These forms significantly improve accuracy of numerical re- 
sults in comparison to direct calculations, like those presented 
in our previous paper [7], see Table [Vl] Nevertheless, with 
256 functions they are about two digits less accurate than the 
most precise results obtained from 9564 Hylleraas basis func- 
tions. These Hyleraas results are slightly more accurate than 
those in |19] due to better optimization of the nonrelativistic 
wave function. 



As we have noticed, for the nonrelativistic energy one needs 
approximately six times smaller basis set of exponentially cor- 
related functions as compared to Hylleraas functions to obtain 
similar accuracy, and the same is confirmed for the mean val- 
ues of operators. The achieved accuracy is limited only by 
the number of basis functions, which nevertheless should be 
well optimized. In practice it demands more computing power 
than we used up to now and a parallel version of the algo- 
rithm would be necessary for optimization of a large number 
of Slater functions. 



I 

to the energy. For these operators we need to use all the dis- 
cussed classes of integrals. g(rii, ri2, na, ri4, ris, ne) integrals 



D. Expectation values 



The basis class of integrals (rii > 0) and the class with 
the one index equal to —1 is sufficient for all mean values 
of operators like those in the Breit-Pauli Hamiltonian, Eq. (|4|i. 
As an example we demonstrate their evaluation for the lithium 
ground state. In Table [VT] we present results for the Dirac 6 
functions using the Drachman formulae iIztIi 



-iEo~V) 



— {En~V) 

Tab 



c ^ 



(19) 



rab 



(20) 



where ^ is a total interaction potential. The similar prescrip- 
tions can be used for operator 



Y^p^) = A{{E, - Vf) ^(Vg^|V^<^). (21) 

a b>c 



In Table IVIII we present numerical values for typical opera- 
tors in higher order perturbation theory i.e. ma^ correction 
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TABLE VII: Expectation value for operators for operators involving extended class of integrals. 



iV 


■s-^ -2 -2 

i. ^ tJ " O 


~2_-2 

i. — tCL^U LtO 


y.-2j,-2 

ii — /ci^fa^'C 0,6 be 


/ ,r,-^U a. ' '^ab'^n nh 

i. id^U oo 


2—12 


1 


204.916 879 


289.283 997 


1.177470 80 


45.945 338 1 


842.720 739 


2 


211.038 011 


297.213 636 


1.112 887 76 


48.080 149 7 


891.514 656 


4 


202.558 151 


284.360 921 


1.07190415 


45.386 279 3 


840.040 129 


8 


202.665 783 


283.643 499 


1.066 132 65 


45.443 462 5 


844.843 155 


16 


202.840 737 


282.053 625 


1.06413154 


45.240 730 3 


853.983 531 


32 


202.852 464 


282.008 017 


1.064188 66 


45.236 709 9 


854.271406 


64 


202.862 382 


282.005 597 


1.064179 55 


45.236 342 7 


854.374 209 


128 


202.883 252 


282.002 183 


1.06417817 


45.236 047 4 


854.554 104 



with two parameters equal to —1 are obtained with double 
integration with 30 x 30 points and all the presented digits 



J 



in Table IVIII are accurate for the corresponding approximate 
wave function. 



r 



III. SUMMARY 

Our primary motivation for developing explicitly correlated 
exponential basis set is the efficient representation of the wave 
function in a small number of basis functions. We applied it 
for the accurate numerical calculation of expectation values 
of some operators corresponding to higher order relativistic 
and QED effects. They involve integrals with quadratic in- 
verse powers of at least two interparticle distances, which are 
the most difficult in the evaluation. Using this compact and 



very flexible correlated exponential basis set, we are aiming 
to determine m and m a' effects in the hyperfine and fine 
structure of lithium-like systems, which have not been inves- 
tigated so far. 
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Appendix A: Weights and nodes of generalized Gaussian 
quadrature with logarithmic end-point singularity 

We present here a set of 30 nodes and corresponding 
weights for generalized Gaussian quadrature with logarithmic 
end-point singularity [26], which was obtained using the al- 
gorithm described in Appendix A of Ref. [19.1 . We used these 
quadrature points for numerical integration of extended inte- 
grals in Tables HVllVll andlVlIl 
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TABLE Vin: Weights and nodes of generalized Gaussian quadrature with logarithmic end-point singularity 



N 


nodes 






weights 




1 


7.323 797 44272605707927651215334608 10^ 


■ a 


2, 


,79892154309547416710987828334736 10" 


-5 


2 


1.100 447 00457774879623368108247943 10^ 


■4 


2, 


.17365526502541548589108155033626 10" 


-4 


3 


5.469 183 26183967432457719790539051 10~ 


■4 


7, 


.20703586534389568501488362593518 10" 


-4 


4 


1.701 857 51910164118701225273639596 10" 


■3 


1, 


.67446096505498972978224436308689 10" 


-3 


5 


4.083 863 60971437462932298732108869 10" 


■3 


3. 


.19128240641146524350794278296664 10" 


-3 




8.300 041 17688233905931162736411746 10" 


■3 


5. 


.35378831352933564469099803155786 10" 


-3 


7 


1.502 297 81560799959114892646895942 10" 


■2 


8. 


.20962136808195567645775658689256 10" 


'3 




2.495 392 36157545503008513656313159 10" 


■2 


1, 


.17680292130848124961801044296255 10" 


-2 


9 


3.878 338 61710629356474590553368630 10" 


■2 


1, 


,59981435048914024332333704241883 10" 


-2 


10 


5.715 089 84811763898261154101387036 10" 


■2 


2, 


,08290410936242947283573768810863 10" 


-2 


11 


8.060 574 14726551690255082886279043 10" 


■2 


2, 


.61515976613093133103470957985343 10" 


-2 


12 


1.095 703 94234517942760556318344690 10" 


■1 


3 


,18220682694563815827694187138199 10" 


-2 


13 


1.443 083 73001500812672361144038109 10" 


-1 


3 


,76672559998067174250314484868819 10" 


-2 


14 


1.848 979 49427531958876501649047717 10" 


■1 


4, 


,34910622632233969988331688489985 10" 


-2 


15 


2.312 130 01548255181588787651651928 10" 


■1 


4, 


,90821532284030403760822595772833 10" 


-2 


16 


2.829 119 60197207748156616199457019 10" 


■1 


5. 


,42224286612626157484616606008227 10" 


-2 


17 


3.394 354 81190852660583946762696458 10" 


■1 


5. 


,86959442909769135605447275253662 10" 


-2 


18 


4.000 131 13109450190847537317246366 10" 


■1 


6. 


,22979181149364796169604799658096 10" 


-2 


19 


4.636 788 56939306293165441422602120 10" 


■1 


6. 


,48434457072330546001338183346875 10" 


-2 


20 


5.292 951 42741036253816689241764505 10" 


■1 


6. 


,61755598512543778100952715354295 10" 


-2 


21 


5.955 843 95325645009469243558808251 10" 


■1 


6. 


,61722952772004782470943263263082 10" 


-2 


22 


6.611 670 40396915198682135573185723 10" 


■1 


6, 


,47524589229209230876914855177978 10" 


-2 


23 


7.246 045 28176032770586180257340226 10" 




6, 


,1879858349954637530119412742859710" 




24 


7.844 457 34734941543894377419518002 10" 


■1 


5, 


.75658036634420432925717153926242 10" 


-2 


25 


8.392 749 51478663018317278487334798 10" 


■1 


5, 


.18697691924656775668123741331021 10" 


-2 


26 


8.877595 97617343115537224494301863 10" 


■1 


4, 


.48981784955672218807632566389701 10" 


-2 


27 


9.286 957 95963227580314671288014994 10" 


-1 


3, 


.68013625003694933496657925550436 10" 


-2 


28 


9.610 500 59187409711054297148665964 10" 


-1 


2, 


.77688703774135920873618253138605 10" 


-2 


29 


9.839 957 03521288981120560285698734 10" 


-1 


1, 


.80238737841607431150476240903211 10" 


-2 


30 


9.969 459 58679763051044061968242492 10" 


■1 


7, 


.82767019549675700264134910161448 10" 


-3 



